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The density-matrix renormalization group is employed to investigate a harmonically-trapped im- 
balanced Fermi condensate based on a one-dimensional attractive Hubbard model. The obtained 
density profile shows a flattened population diff'erence of spin-up and spin-down components at the 
center of the trap, and exhibits phase separation between the condensate and unpaired majority 
atoms for a certain range of the interaction and population imbalance P. The two-particle density 
matrix reveals that the sign of the order parameter changes periodically, demonstrating the real- 
ization of the Fulde-Ferrell-Larkin-Ovchinnikov phase. The minority spin atoms contribute to the 
quasi-condensate up to at least P ~ 0.8. Possible experimental situations to test our predictions 
are discussed. 

PACS numbers: 03.75.Hh, 71.10.Pm, 71.10.Fd 



CD 



O 

o 



> 

00 

q 

00 

o 
l> 
o 



- '—I 

X 
J3 



Over the past few years, two major breakthroughs have 
been achieved in fermionic superfluids of tenuous atomic 
vapor: BEC-BCS crossover P, 0] and imbalanced super- 
fluidity d, 3- These subjects have been studied not just 
in atomic physics but also in diverse subfields of physics 
such as condensed matter physics d, 0, H and nu- 
clear physics El, [13] ■ Two major issues in imbal- 
anced superfluidity are whether superfluidity disappears 
at a particular value (the Chandrasekhar-Clogston (CC) 
limit [isj ) of population imbalance, and in what parame- 
ter regime the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
phase [3, emerges. The observation of the CC limit 
is currently under controversy, while the FFLO phase 
remains elusive flil despite extensive theory literature 



H H 0, EUii, 0, [11, 111, [i^ 

In this Letter, we address these issues by appl ying the 
density-matrix renormalization group (DMRG) [U 
to a system of harmonically trapped fermions. We con- 
sider a trapped one-dimensional (ID) system, which can 
exhibit BEG due to the cut-off of infrared divergence. 
We find that in the ground state of an imbalanced Fermi 
system the pairing order parameter undergoes a pe- 
riodic sign change — the hallmark of the FFLO state 
29l . 30l . 3l| . We also show that, as the imbalance becomes 
greater, the minority-spin atoms continues to contribute 
to the quasi-condensed state, implying non-existence of 
the CC hmit. 

We consider a system of spin- 1/2 fermions undergoing 
short-ranged interaction and confined in a ID harmonic 
potential whose characteristic size is We take I as the 
unit of length and discretize the system by introducing a 
lattice of L sites with the lattice constant given by d = 
21/ L. The transfer amplitude between neighboring sites, 
t — li^ /2m(P (m is the mass of the atom), reproduces 
the energy dispersion of the free space in the limit of 
small filling factor {L oo). The contact interaction 



giB^iz^ — zi) is approximated by introducing an on-site 
interaction with coupling constant U — gin/d between 
atoms in different spin states. 

We use DMRG to calculate the on-site pair correla- 
tion function and the two-body reduced density matrix 
(2BDM) for the L-site Hubbard model with a harmonic 
on-site potential. DMRG allows an efficient, numerically 
exact treatment of many-body problems in ID by itera- 
tive truncation of the Hilbert space. We retain 200 trun- 
cated states per DMRG block with the maximum trun- 
cation error of 10~^. 

The Hamiltonian of our system is given by 
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where c^.o- annihilates an atom at site i in spin state cr(=t 



a), 



fii^l -\- fii^i, and V = 4A/L'^ 



is determined from the depth A of the potential. We 
note that A/t oc and U/t oc L~^. We calculate 
the ground state of Hamiltonian ([T|) within the particle- 
number sector of N-^ and atoms in spin states f and |, 
respectively. The imbalance parameter is defined as P = 
(iV| -Ni)/N, where N = -\-Ni . The Fermi momentum 
at the trap center is calculated from the averaged density 
Tier = {hi^cr)L/2l as /efct = riaTr, where the overbar denotes 
the average over O.IL ~ {).2L neighboring sites. 

The s-wave scattering length of atoms ai^ in a ID 
trap with radial width a±_ = ^Jh/fiuj^, where /i is half 
the atom mass and lo_\_ is the radial trapping frequency, 
is modified from the free-space scattering length osd 
[s^ ]: in the low-energy limit of incident atoms, aro — 
-{a\/2ay.Yi) X (1 - Ca3D/a_L) with C ~ 1.46. Thus the 
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(Color online) On-site pair correlation 
(left column) and density distributions 
(right column) of spin-up and spin-down atoms together 
with their difference plotted against the z coordinate in a 
harmonic potential for {A/t,U/t) = (6400/L^ -800/L) and 
{Nt,Ni) = (a) (20,20), (b) (22,18) and (c) (28,12) with an 
L — 200-site lattice. In (a) the density of spin-up atoms, 
and in (b) and (c) the difference for the noninteracting case 
((7 = 0) are also plotted. 



ID effective interaction is described by U(z) — (7id<5(2;), 
where gm = -h^/tiaiu- We take A/t = 6400/L^, 
which, in the case of ^"K atoms with a± = 86 nm and 
uj, = 2ttx 256 Hz [SSj], gives I = 6.28 fim; U/t = -800/L, 
for example, corresponds to aiD ~ 62.8 nm. 

We first examine the on-site pair correlation function 
defined by 



O. 



on — site 
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where ji/'o) is the ground state of the system. The left 
column of Fig. [1] displays Oon-sito(2j, zj) for three values 
of P with L = 200 and U/t = -4 [(fepaio)"^ = 1.82 
for n| = = 17.5]. Figure [TJa) exhibits the case with 
P = 0, where Oon-sito(^i, Zj) shows a slow decay without 
changing the sign, and then drops precipitously towards 
zero at \zi\ ~ 0.75. The sharp boundary refiects cohesion 
of the system due to attractive interaction. Figure \V[h) 
shows the case with P = 0.1, where the pair correlation 
function changes its sign in real space, and the amplitude 
of the oscillation vanishes rather sharply at \zi\ ^ 0.7. 
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FIG. 2; (Color online) (a) Dependence of the difference in 
density between spin-up and spin-down atoms, — ni, on z 
calculated for A'^ = 40 atoms with imbalance parameter P = 
0.1, L = 64, 80, 128, 200 and {A/t, U/t) = (6400/L^ -800/L). 
n| — ni for P = 0.2 with L — 200 is also plotted, (b) 
Wave vector of oscillations in the pair correlation and that 
in n| — ni obtained for = 40 atoms for various values 
of the imbalance parameter with L = 200 and [A/t, U/t) = 
(0.16,-4). (c) Phase diagram for A/t = 6400/L^ Hi] LO 
denotes the Larkin-Ovchinnikov state. Here, fcp is defined as 
(fcpi +A;f|)/2. Phase separation occurs regardless of the value 
of (fcFaiD)^^ when P > 0.21 (dotted line). 



Figure [T] (c) shows the case with P — 0.4, where the 
sign change of the pair correlation function occurs with 
a much shorter period because of a larger separation in 
the Fermi wave numbers of the up and down spin compo- 
nents. The peaks of Oon-sitdzi, Zj) align along straight 



lines of constant 



implying that the order param- 



eter oscillates periodically in space. In Figs. [T] (b) and 
(c), the order parameter varies sinusoidally with no bro- 
ken time-reversal symmetry, indicating that the realized 
states are in the Larkin-Ovchinnikov (LO) phase rather 
than the Fulde-Ferrell phase. 

The right column of Fig. [T] shows the density distri- 
butions of the spin-up and spin-down atoms together 
with their difference. For P = 0, the total number of 
atoms peaks at the center of the harmonic trap with small 
Friedel oscillation due to the presence of the trapping po- 
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tential. We note that the pair correlation extends over 
the whole region where atoms are present. In Fig. [T] (b) 
with P — 0.1, the difference in site occupation numbers is 
almost flat near the trap center. The population plum- 
mets again at the locations where the pair correlation 
function vanishes. In Fig. [1] (c) with P = 0.4, we see 
that the tails of the total population at both ends only 
comprise the majority atoms in the | t) state. The pop- 
ulation difference first increases rapidly as we go from 
an edge toward the center, and then becomes almost flat 
near the trap center. The oscillating pair correlation in 
Fig. [1] (c) disappears at \zi\ ~ 0.6 at which the number 
of the minority population vanishes. This fact indicates 
that the two phases — the FFLO condensate with almost 
constant density of excess majority atoms and the nor- 
mal, spin-polarized state — phase-separate. 

Figure [5] (a) plots the difference of density distribu- 
tions of spin- up and spin-down atoms for P = 0.1 with 
varying L, and for P = 0.2 with L = 200. We find 
that the density difference shows a rapid convergence, 
which indicates that L — 200 is close to the continuum 
limit, and that the difference oscillates around a slowly- 
varying parabolic curve around the center of the trap. 
Near the trap center, the oscillation of the density dif- 
ference is incommensurate with the lattice and almost 
independent of the lattice constant for L > 80. By a 
nonlinear fit we determine the wave vector of the oscil- 
lation as well as the density difference at z = (the 
peak of the parabola). Figure [2] (b) plots the wave vec- 
tor of the oscillations against the wave-vector difference 
Afcp = fcpT — k-pi for the pair correlation function and 
for the density difference — at the trap center. The 
linear relation q = Afcp, which holds for the former case, 
is consistent with the LO phase [2^, 3^, 31 1 . The relation 
q = 2 Afcp, which holds for the latter case, is expected for 
the FFLO states [sl], and this is confirmed in Fig.[2](b). 

The sudden rise of the minority population for P > 
does not appear in the noninteracting system and show 
striking resemblance to the density profiles observed in 
the Rice experiment [36|. Such cohesion in the density 
distribution implies that the minority atoms are drawn 
into the inner core due to pair correlation, while un- 
paired atoms in the majority state are pushed outside 
the core. While similar density profiles have been repro- 
duced in 3D simulations with phenomenological surface 
tension 37|, |38| , it is interesting to observe that such a 
steep rise also occurs in the ID system by treating many- 
body effects rigorously. We identify the point of phase 
separation with the onset of the two shoulder peaks in 
the density difference. Figure [2] (c) shows the phase di- 
agram for A/t — 6400/L^ obtained by identifying the 
onset for 36 < < 41. As the population imbalance P 
becomes greater, stronger on-site attractive interaction 
or shorter aiD is needed to eliminate phase separation, 
and for P > 0.21, the shoulder peaks are always seen 
regardless of the strength of the on-site attractive inter- 



CO 



(a) 



P=0.00 
P=0.10 
P=0.20 
P=0.50 
P=0.80 
P=0.90 




10 20 30 
Eigenstate number 




0.2 0.4 0.6 0.8 1 
Imbalance parameter P 



-9 20 



o 

< 10 
_^ 

re 5 



.1, 







^^^^^•-...^ minonty atoms 




sum of large eigenvalues*"'^':^:!^ 








U80 {A/t, U/f)=(1, -10) /V=40 





0.2 0.4 0.6 0.8 1 

Imbalance parameter P 

FIG. 3: (Color online) (a) Eigenvalue distribution for the two- 
particle density matrices calculated for various values of the 
imbalance parameter P for {A/t,U/t) = (6400/L^ -800/L) 
and A'^ = 20 on a L = 80-site lattice, (b) Maximum eigenvalue 

^Max f^o'" different total numbers of atoms A'^ — 20, 30, 40. (c) 
The sum of eigenvalues up to the kink, and the number of the 
minority spin atoms plotted against P. 



action. 

Finally, we study the dependence of the condensate 
fraction on imbalance parameter P by calculating the 
eigenvalue distribution of the two-body density matrix 
(2BDM) [sil. The matrix elements of 2BDM in the site 
representation are given by 



(2) 
Pii' 



.TCjuI^o)- 



(3) 



Note that p^-'-j equals Oon-sitc(«, j)- The sum of the 
eigenvalues equals N-^N^. For the noninteracting case, 
there is an iV|iV|-fold degenerate eigenvalue of unity, and 
the rest of the eigenvalues are all zero for the ground 
state, because each single-particle energy level is either 
occupied or unoccupied. 

Figure [3] (a) shows the distribution of the eigenvalues 
of 2BDM for various values of the imbalance parame- 
ter P. We observe that the first few eigenvalues stand 
out from the rest. This can be interpreted as the atoms 
forming a quasi-condensate, where pairs of atoms con- 
dense into more than one eigenstate. Figure [3] (b) shows 
how the maximum eigenvalue decreases with increasing 

(2) 

P. It also shows that ej^^^^ is rather insensitive to the 
value of N. This fact suggests that the lowest-lying pair- 
ing state is already maximally occupied for a small value 
of N and that additional atoms contribute to other pair- 
ing states with oscillating sign changes compatible with 
the LO state. Figure [3] (c) shows the sum of those large 
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eigenvalues, together with the number of the minority 
atoms. We note that the sum is almost equal to the 
number of minority-spin atoms, indicating that almost 
all of the latter contributes to the quasi-condensate. This 
phenomenon should be interpreted as being due to the 
abundance of majority-spin atoms which can pair with 
minority-spin atoms no matter where the latter reside. 
This behavior and the slow decay of the maximum eigen- 
value of 2BDM suggests the robustness of condensate 
against population imbalance. 

The eigenvector of the 2BDM with the largest 

f 2) 

eigenvalue e^ax' which describes the state in which the 
largest number of Cooper pairs participate, is contributed 
mostly from the i — j components. This reflects the fact 
that the two atoms on the same site are most effectively 
paired, as expected from Hamiltonian ([T]). The sign of 
changes as a function of i when P > 0, showing that 
the sign change in the order parameter originates from 
the FFLO nature of the condensate rather than from the 
Andreev scattering of particles at the normal-superfluid 
boundary 26]. 

In conclusion, we have employed the density-matrix 
renormalization group to investigate fermionic conden- 
sation of a ID trapped gas with population imbalance 
between the two spin states, and have shown that the 
ground state exhibits the FFLO-like feature over a wide 
range of imbalance parameter P. We have also found 
a critical value of P beyond which the LO condensate 
and spin-polarized Fermi gas phase-separate at T = no 
matter how strong the interaction is. 
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